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Abstract 

We consider the scattering of entangled two-photon states from collections of small particles. 
We also study the related Mie problem of scattering from a sphere. In both cases, we calculate the 
entropy of entanglement and investigate the influence of the entanglement of the incident field on 
the entanglement of the scattered field. 
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The propagation of quantum states of light in complex media such as the atmosphere, 
colloidal suspensions and biological tissue is a topic of fundamental interest and considerable 
applied importance. New physical phenomena, including the transmission of quantum states 
through random media [1~4]; the observation of spatial correlations in multiply-scattered 
squeezed light [5, 6]; and the measurement of two-photon speckle patterns [7-9] have been 
described. Applications to spectroscopy [10], imaging [11-22] and communications [23-27] 
have also been reported. To understand the underlying physical principles, it is often useful 
to consider relatively simple model systems. A step in this direction was taken in [28], 
where the propagation of two-photon entangled states in random media was studied within 
the framework of radiative transport theory. In this Letter, we investigate the scattering of 
entangled two-photon states in an even simpler setting, namely from deterministic media. 
We study in detail the problem of scattering from collections of small particles, as well as 
the related Mie problem of scattering from a sphere. In particular, we analyze the extent 
to which scattering can alter the quantum correlations of the optical held and calculate the 
resulting entropy of entanglement. 

We begin by considering the propagation of a quantized held in a material medium with 
dielectric permittivity e. For simplicity, we work with the scalar theory of the electromag¬ 
netic held. The electric-held operator E may be decomposed into positive and negative 
frequency components: E = E^ + E~. The positive frequency component E^ obeys the 
wave equation [29, 30] 


V^E+ = 


e r 


d‘^E+ 


BE 


( 1 ) 


The negative frequency component E~ is dehned by E~ = [77"'']^. Here the medium is taken 
to be nonabsorbing, so that e is purely real and positive. Now, let |'^) be a two-photon state 
and dehne the second-order coherence function T*^^) as the normally ordered expectation of 
held operators: 


r(2)(r,f;r ,f') = (^1(r,t).F (r',t') 

xE^{r',t')E^{Y,t)\'ip) . 


( 2 ) 


We note that T^^^ is proportional to the probability of detecting one photon at r at time 
t and a second photon at r' at time t’. It can be measured in a Hanbury-Brown-Twiss 
interferometer [31]. It can be seen that T^^^ factorizes [32] as 

T^^\r,t;r,t') = \A{r,t;v ,t')\^ . (3) 
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FIG. 1. Illustrating the notation. 


Here |0) is the vacuum state and the two-photon amplitude A is dehned by 

A{r, t] r, t') = (0| E^{y, t).E+(r', t') \%l)) . 


( 4 ) 


Evidently, A satishes the pair of wave equations 


VlA 

Vl,A 


e{r) d‘^A 

’ 

e(r') d^A 
dt'"^ ’ 


( 5 ) 

( 6 ) 


which follow from the fact that obeys the wave equation [33]. We will hud it convenient 
to introduce the Fourier transform of the two-photon amplitude A, which is dehned by 


A{i 




Eq. (5) then becomes 


( 7 ) 


VlA + kh{r)A = 0 , 

( 8 ) 

Vl,A + E^e{E)A = 0 , 

( 9 ) 


where k = u/c and k' = ca'/c. 

We now develop the scattering theory for the two-photon amplitude. To proceed, we 
consider the Helmholtz equation 


-I- /c^e(r)n = 0 . 


( 10 ) 


The held u is taken to consist of incident and scattered parts, which we write as the sum 
u = Ui + Us- The incident held Ui is the held that would exist in the absence of the scatterer. 
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The scattered field Us is given by [34] 

M^(r) = j d^ri£r2G{T,Ti)T{ri,r2)ui{T2) , (11) 

where the Green’s function G is of the form 

ik\r-r'\ 

= ( 12 ) 

The T-matrix obeys the integral equation 

T(r,r') = k‘^ri{r)6{r — r') + k^ri{r) j (i^r"G(r, r")T(r", r') , (13) 

where the susceptiblity rj is defined by the relation e = 1 + Anrj. Applying (11) to each 
of the arguments of A, we find that the Fourier transformed two-photon amplitude of the 
scattered field, denoted is given by 

A5(r,r') = J d^ri£r2d^r[£r2G{r,ri)T{ri,r2) 

xG'(r',r'jT'(r),r2)Ai(r2,r2) , (14) 

where Aj is the Fourier-transformed two-photon amplitude of the incident field. In addition, 
G' and T' denote the Green’s function and T-matrix at the wavenumber k'. If As(r,r') 
factorizes into a product of two functions which depend upon r and r' separately, we will 
say that the two-photon state \%Ij) is not entangled. In contrast, an entangled state is not 
separable. We can now state our first result. It follows directly from (14) that if Aj is 
separable then A^ is separable. That is, if A^ is entangled then Aj is entangled, which 
means that entanglement cannot be created by scattering an unentangled incident state. 

We now consider the far-field limit of the two-photon amplitude A 5 (r,r'), where the 
points of observation r and r' he on a sphere of radius R in the far-zone of the scatterer, 
as shown in Figure 1. In doing so, we make use of the asymptotic behavior of the Green’s 
function G(r, r') for r 3> r': 

pikr 

G(r,r')-. (15) 

r 

We also use the plane-wave expansion for the two-photon amplitude Aj, which is of the form 

Aj(r,r') = J dkdk'A(k,k')G(^*^'"+^'*^'-'') . (16) 

Here Aj is a suitable coefficient, which is expressible as a vacuum to two-photon state 
transition amplitude. Thus (14) becomes 

A(k, k') = J dkidks (k| T |ki) (k'l T' Iks) Aj(ki, ks) , (17) 
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where k = /cr, k' = k'i' and A = exp [—i{k + k')R] Ag, with |r| = |r'| = R. See Figure 1 
where the geometry of the problem is illustrated. We note that the above T-matrices are 
on-shell. The momentum-space T-matrix elements are dehned by 

(k|T|k') = J r) , (18) 

where |k| = |k'| = k. Given a scattering medium characterized by its T-matrix, (17) 
predicts the two-photon amplitude of the far-zone scattered held in terms of the two-photon 
amplitude of the incident held. As may be expected, if Ai is separable then A is separable, 
consistent with (14). If ^j(ki,k 2 ) = (5(ki — £ 2 ), which corresponds to a fully entangled 
two-photon state, then (17) becomes 

^(k, k') = J dk" (k| T |k") (k'l r |k") , (19) 

where |k| = k and |k'| = k'. 

We now compute A for several diherent scattering systems. We begin with a small 
spherical scatter of radius a, where ka 1. The T-matrix is then given by 

(k|T|k') =t(A;)e*(‘^-‘^')'"« , (20) 

where fq is the position of the scatterer and t{k) is dehned in the supplementary material. 
Making use of (19) and (20) we obtain 

.4(k,k') (!(«: +/t')r„|) . 

( 21 ) 

We see at once that A is separable and thus the scattered held is unentangled, even when 
the incident held is entangled. 

Next, we consider a collection of identical small scatterers. The T-matrix is given by 

(k|T|k') = 5^U(A:)G(*^'^“-i^'-^*>) , (22) 

a,b 

where {ra} are the positions of the scatterers and tab is dehned in the supplementary material. 
Using (19), we hnd that A is given by 

M(k,k') = 

a,b a',b' 

xsinc(|/crb-|-/cTfoil) . (23) 
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FIG. 2. Two-photon amplitude .A(k, k') of a spherical scatterer of radius a as a function of the 
angle 8 between k and kh The index of refraction of the sphere is n = 1.5. 

We note that in general A is nonseparable; thus the scattered held is entangled. 

Finally, we consider a homogeneous spherical scatterer of radius a centered at the origin 
with index of refraction n. The T-matrix is of the form 

(k| T |k') = 5^(2/ + l)Ai{k)Piik ■ k') , (24) 

i 

where the Mie coefficient Ai is dehned as [35] 

1 ji{nka)j[{ka) - nji{ka)j[{nka) 

-7 /-|\ /-iW • 

nh\ \ka)j[{nka) — h\ {ka)ji{nka) 

Using the identity 

j dk!'Pi{k ■ k")PKk' ■ k") = ^^^5wPi{k ■ k') (26) 

to carry out the integral in (19), we hnd that A is given by 

^(k, k') = dvr 5^(2/ + l)Ai{k)Ai{k')Pi{k • k') . (27) 

i 

In Figure 2 the quantity ^(k, k') is plotted as a function of the angle between k and k'. 

We now turn to the computation of the entanglement entropy for the above systems. 
Following [36], we consider the singular value decomposition (also known as the Schmidt 
decomposition) of the two-photon amplitude, viewed as an operator with kernel ^(k, k'). 
We hnd that ^(k, k') can be decomposed into a superposition of separable terms of the form 

‘^(k,k') = ^anUn{k)vl{k!) , (28) 
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where each term can be interpreted as not entangled. Here the singular values an are real¬ 
valued and the singular functions obey 


J {A*A) (k,k')t;„(k')dk' = a^n„(k) , 

(29) 

1 {AA*) (k, k')M„(k')(ik' = a^Mn(k) , 

(30) 


where A* is the adjoint of the operator A. A measure of the degree of entanglement is the 
entropy S, which is dehned by 

S = -^an\ogan. (31) 

n 

We note that the larger the value of S', the greater the degree of entanglement. 

To illustrate the above results, we calculate the entanglement entropy of a spherical 
scatterer. Using (27) and the identity (26), we find that 

(A'A) (k, k') = (to)* 5^ (Ai(k)A,{k')f K,.,.(k)y);(k') . (32) 

l,m 

We immediately see that the singular functions and singular values are given by 

Ulm{k) = Vlm{k) = Um(k) , (33) 

(^Im = cri = {^7rf\Ai{k)Ai{k')\ . (34) 

Thus the entropy is given by the formula 

S =-'^{21 + 1) ai log ai . (35) 

i 

In Fig. 3 we plot the entropy as a function of the radius of the sphere. We see that in the 
limit where the radius tends to zero, the entropy vanishes, consistent with the separability 
of the two-photon amplitude for the case of a point scatterer in (20). We note that the 
entropy is oscillatory and increasing, but not monotonically. Thus, large spheres generally 
have greater entropies than small spheres. We also note that the oscillations are related to 
the presence of scattering resonances in the Mie coefficients Ai. 

Next, we consider a system of point-scatterers. Using (23), the identity 

= dTT ^ ii{xy)Yi^{^)Y^^{y) (36) 

l,m 

and the orthogonality of the spherical harmonics, we see that A*A can be written in the 
form 
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FIG. 3. Entropy of entanglement of a spherical scatterer as a function of radius a. The index of 
refraction of the sphere is n = 1.5. 


where 


{A*A) (k, k') 


l'm'^ 


l,m I',m' 


(37) 


A 


Im 

I'm' 


E ^lm* ^l"m" 
^l"m"^l'm' 

l",m" 


(38) 


The coefficients Cf^' contain the information on the positions of the scatterers and are 
defined by 


G/,™, = (dvr)^ ^i'+''Gb(fc)W(A;')sinc {\kra + fc'r„/|) 

a, a' 

b, b' 


(39) 


To construct the singular value decomposition of A, we expand the singular functions Un 
(which satisfy (29)) into spherical harmonics of the form 


Wrx(k) = ’ 

l,m 


(40) 


where the coefficients are to be determined. Making use of (37) and the orthogonality 
of the spherical harmonics, we find that the can be obtained from the solution to the 
eigenproblem 


E 

I' ,m' 


A 


l'm'(n) _ _2,,(7i) 


Im ^I'm' 


= <^rAlm 


(41) 



Once the above eigenproblem has been solved, the entropy is computed from (31). In Fig. 4 
we plot the entropy as a function of the distance between a pair of point scatterers. We see 
that the entropy decreases as the separation between the scatterers increases. In the limit 
where the scatterers are far apart (noninteracting), the entropy vanishes, consistent with 
our results for the case of a single point scatterer. 

We close with a few remarks, (i) There is a well-known duality between partially coherent 
and partially entangled light [37]. We note that our results are analogous to the change in 
coherence that can occur with scattering [31]. (ii) It has been reported that entanglement can 
be induced by multiple scattering in random media [4]. The opposite conclusion was argued 
in [28]. Since we consider deterministic systems, the results of this Letter are potentially 
signihcant because they remove from consideration the role of randomness in modifying the 
entanglement of the incident held, (hi) Although in our model the electromagnetic held 
is quantized, the interaction of the held with the scattering medium is treated classically. 
It would be of interest to extend our results to the case in which the medium consists of 
a collection of two-level atoms. In this manner, it should be possible to understand the 
transfer of entanglement from the held to the medium [38]. 

In conclusion, we have studied the scattering of entangled two-photon states from non¬ 
absorbing material media. In the setting of simple model systems, we have calculated the 
entropy of entanglement and have characterized the inhuence of the entanglement of the in¬ 
cident held on the entanglement of the scattered held. In future work, we plan to investigate 
whether the observed resonances in the entropy are statistically stable in random media. 



kd 


FIG. 4. The entropy of entanglement of a pair of scatterers separated by a distance d. The radii 
of the scatterers is ka = 0.2 and their index of refraction is n = 1.5. 


9 



The authors are grateful to Paul Berman for valuable discussions. This work was sup¬ 
ported in part by the NSF Center for Photonic and Multiscale Nanomaterials under the 
grant DMR-1120923. 


[1] P. Lodahl, A.P. Mosk and A. Lagendijk, Phys. Rev. Lett. 95, 173901 (2005) 

[2] P. Lodahl and A. Lagendijk, Phys. Rev. Lett. 94, 153905 (2005) 

[3] M. Cande, A. Goetschy and S. E. Skipetrov, Europhys. Lett. 107, 54004 (2014) 

[4] J. R. Ott, N. A. Mortensen and P. Lodahl, Phys . Rev . Lett. 105, 090501 (2010) 

[5] S. Smolka, A. Huck, U. L. Andersen, A. Lagendijk and P. Lodahl, Phys. Rev. Lett. 102, 
193901 (2009) 

[6] S. Smolka, J. R. Ott, A. H. Ulrik, L. Andersen and P. Lodahl, Phys. Rev. A 86, 033814 (2012) 

[7] W. H. Peelers, J. J. D. Moerman and M. P. van Exter Phys. Rev. Lett. 104, 173601 (2010) 

[8] H. D. Pires, J. Woudenberg and M. P. van Exter, Phys. Rev. A 85, 033807 (2012) 

[9] C. W. J. Beenakker, J. W. F. Venderbos and M. P. van Exter, Phys. Rev. Lett. 102, 193601 
(2009) 

[10] S. E. Skipetrov, Phys. Rev. A 75, 053808 (2007) 

[11] D. N. Klyshko, Zh. Eksp. Teor. Fiz. 94, 82 (1988) [Sov. Phys. JETP 67, 1131 (1988)] 

[12] D. V. Strekalov, A. V. Sergienko, D. N. Klyshko and Y. H. Shih, Phys. Rev. Lett. 74, 3600 
(1995) 

[13] A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko and M. C. Teich, Phys. Rev. Lett. 87, 
123602 (2001) 

[14] A. F. Abouraddy, P. R. Stone, A. V. Sergienko, B. E. A. Saleh and M. C. Teich, Phys. Rev. 
Lett. 93, 213903 (2004) 

[15] A. Gatti, E. Brambilla, M. Bache and L. A. Lugiato, Phys. Rev. Lett. 93, 093602 (2004) 

[16] G. Scarcelli, A. Valencia and Y.H. Shih, Europhys. Lett. 68, 618 (2004) 

[17] G. Scarcelli, V. Berardi and Y.H. Shih, Phys. Rev. Lett. 96, 063602 (2006) 

[18] B. 1. Erkmen and J. H. Shapiro, Phys. Rev. A 78, 023835 (2008) 

[19] M. D’Angelo, A. Valencia, M.H. Rubin and Y.H. Shih, Phys. Rev. A 72, 013810 (2005) 

[20] J. C. Schotland, Opt. Lett. 35, 3309 (2010) 

[21] M. B. Nasr, B. E. A. Saleh, A. V. Sergienko and M. C. Teich, Phys. Rev. Lett. 91, 083601 


10 



(2003) 

[22] M. C. Teich, B. E. A. Saleh, F. N. C. Wong and J. H. Shapiro, Quant. Inf. Process. 11, 903 

( 2012 ) 

[23] Zhen-Sheng Yuan, Xiao-Hui Bao, Chao-Yang Lu, Jun Zhang, Cheng-Zhi Peng and Jian-Wei 
Pan, Phys. Rep. 497, 1 (2010) 

[24] J. Tworzydlo and C. W. J. Beenakker, Phys. Rev. Lett. 89, 043902 (2002) 

[25] A. L. Moustakas, H. U. Baranger, L. Balents, A. M. Sengupta and S. H. Simon, Science 287, 
287 (2000) 

[26] S. E. Skipetrov, Phys. Rev. E 67, 036621 (2003) 

[27] J. H. Shapiro, IEEE J. Selected Topics in Quantum Electronics, 15 (2009) 

[28] V. A. Markel and J. C. Schotland, Phys. Rev. A. 90, 033815 (2014) 

[29] R. J. Glauber and M. Lewinstein, Phys. Rev. A 43, 467 (1991) 

[30] S. Scheel, L. Knoll, D.-G. Welsch and S. M. Barnett, Phys. Rev. A 60, 1590 (1999) 

[31] L. Mandel and E. Wolf, Optical Coherence and Quantum Opties (Cambridge University Press, 
Cambridge, 1995) 

[32] M. H. Rubin, Phys. Rev. A 54, 5349 (1996) 

[33] B. E. A. Saleh, M. C. Teich and A. V. Sergienko, Phys. Rev. Lett. 94, 223601 (2005) 

[34] E. N. Economou, Green’s Functions in Quantum Physics, third edition (Springer Verlag, 

Berlin, 2006) 

[35] W. T. Grandy, Scattering of Waves by Large Spheres (Cambridge University Press, Cambridge, 
2005) 

[36] C. K. Law and J. H. Eberly, Phys. Rev. Lett. 92, 127903 (2004) 

[37] B. E. A. Saleh, A. F. Abouraddy, A. V. Sergienko and M. C. Teich, Phys. Rev. A 62, 043816 

( 2000 ) 

[38] P. R. Berman, Phys. Rev. A 76, 042106 (2007); ibid. 76, 043816 (2007) 


11 



Supplementary Information 


Here we collect some basic results about point scatterers. First, consider a small spherical 
scatterer of radius a, where ka 1. The T-matrix is then given by 

(k|T|k') , (42) 

where t{k) = ak'^ and ro is the position of the scatterer. The renormalized polarizability a 
is dehned by 

^ _(^0_ / . o '! 

1 — (2a) — iaok'^ ’ 

where the polarizability oq = a^(n^ — l)/3, with n the index of refraction. Note that the 
above formula includes radiative corrections to the Lorentz-Lorenz form of the polarizability. 
Next, we consider a collection of identical point scatterers. The susceptibility is of the form 
rjir) = ? 7 o Yla where rjo = {v? — l)/47r, {ra} are the positions of the scatterers and 



The T-matrix is given by 


1 , |r| < a , 
0 , |r| > 0 . 


where 


Here 


(k|T|k') = 5^U(A:)e*(*^'^“-*^''^^) , 

a,b 

tab = aok‘^M~^ . 


tkkab ^ab j 


where 



G(ra,rb) , a b , 

3 

-h ^k , a = b . 

2a 


(44) 

(45) 

(46) 

(47) 

(48) 
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